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Dirac fermions have been studied intensively in condensed matter physics in recent years. 
Many theoretical predictions critically depend on the number of valleys where the Dirac 
fermions are realized. In this work, we report the discovery of a two dimensional system 
with a single valley Dirac cone. We study the transport properties of HgTe quantum wells 
grown at the critical thickness separating between the topologically trivial and the quantum 
spin Hall phases. At high magnetic fields, the quantized Hall plateaus demonstrate the 
presence of a single valley Dirac point in this system. In addition, we clearly observe the 
linear dispersion of the zero mode spin levels. Also the conductivity at the Dirac point and 
its temperature dependence can be understood from single valley Dirac fermion physics. 



I. INTRODUCTION 



In recent years, Dirac fermions have been intensively studied in a number of condensed matter 
systems. In the two dimensional material graphene the low energy spectrum is well described by 
two spin degenerate massless Dirac cones at two inequivalent valleys, giving rise to four mass- 
less Dirac cones in totaU* 2 -. The fabrication of graphene sheets enabled substantial experimental 
progress in this field, and the physics of the Dirac fermions has been investigated extensively-. At 
the same time, many theoretical predictions rely on a single Dirac cone valley, or, at least, weak 
inter- valley scattering 3 . Graphene is not a suitable platform to test these latter predictions because 
of the presence of two valleys and strong inter- valley scattering. In addition, it is presently unclear 
how an energy gap can be reliably generated in single layer graphene, which would be desirable 
for a variety of device applications. 

A HgTe/CdTe quantum well is another system where Dirac fermion physics emerges^. In 
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this case, the Dirac fermions appear only at a single valley, at the T point of the Brillouin zone. 
Furthermore, tuning the thickness d of the HgTe quantum well continuously changes both the 
magnitude and the sign of the Dirac mass. When d is less than a critical thickness d c ~ 6.3 nm, the 
system is in a topologically trivial phase with a full energy gap. On the other hand, when d > d c , 
the quantum spin Hall state is realized, where a full energy gap in the bulk occurs together with 
gapless spin-polarized states at the edge. The experimental discovery of this state 5 provides the 
first example of a time-reversal invariant topological insulator in nature 6 . A topological quantum 
phase transition is predicted to occur when d = d c , where a massless Dirac fermion state is realized 
at a single valley, with both spin orientations 4 . Our paper reports the experimental discovery of 
such a state. 

In a two dimensional system with time reversal symmetry and half integral spin, a minimal 
number of two massless Dirac cones can be present, as can be proven by a simple generalization 
of a similar theorem in one dimension 7 . General results of this type have first been discovered in 
lattice gauge theory, and are known as chiral fermion doubling theorems 8 . In this sense, the HgTe 
quantum well at the critical thickness d = d c realizes this minimal number of two Dirac cones in 
two dimensions 9 . 

II. HGTE QUANTUM WELLS AS HALF-GRAPHENE 

HgTe is a zinc -blende-type semiconductor with an inverted band structure. Unlike conventional 
zinc-blende semiconductors, and due to the very strong spin-orbit coupling in the material, the T 8 
band of HgTe (which derives from chalcogenide p-orbitals), has a higher energy than the T 6 band 
that originates from metallic s-orbitals and usually acts as the conduction band. Consequently, in 
HgTe/(Hg,Cd)Te quantum wells, when the well thickness is large enough, the sub-bands of the 
quantum well are also inverted: T 8 derived heavy hole-like (H) sub-bands have higher energies 
than r 6 -based electron-like (E) sub-bands. The inverted band structure, especially the inversion 
between E 1 and Hi sub-bands (where the suffix is the sub-band number index), leads to the oc- 
currence of the quantum spin Hall effect 4 - 5 , boasting dissipationless edge channel transport at zero 
external magnetic field 10 . When the thickness of the quantum well is decreased, the energies of 
the E sub-bands increase due to quantum confinement, while those of the H sub-bands decrease, 
as shown in Fig[T](a). Eventually, the E\ sub-band gains a higher energy than the H\ sub-band 
and the system has the normal band sequence. The different dependence of E and H sub-bands 
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on well thickness implies that there must exist a critical thickness where the band gap is closed. In 
fact, the crossing point between E 1 and Hi sub-bands, denoted as d c in Fig[T](a), not only corre- 
sponds to the critical point for the quantum phase transition between quantum spin Hall insulator 
and normal insulator 4 but additionally yields a quantum well whose low-energy band structure 
closely mimics a massless Dirac Hamiltonian. 

When d = d c ~ 6.3 nm, the energy dispersion of the E\ and Hi sub-bands, which can be 
calculated from an 8-band Kane model, is found to linearly depend on the momentum k near the 
T point of the Brillouin zone, as shown in Fig Q] (b) and (c). Near the T point, using the states 
\Ei, |), l-ffi, |), \Ex,~) and \H X , — |) as a basis, one can write an effective Hamiltonian for the 
Ei and Hi sub-bands, as follows 4 



The two components of the Pauli matrices a denote the Ei and Hi sub-bands, while the two 
diagonal blocks H D (k) and H* D {— k) of H e ff represent spin-up and spin-down states, related to 
each other by time reversal symmetry. At the critical thickness, the relativistic mass M. in (O 
equals to zero. If we then only keep the terms up to linear order in k for each spin, Ho(k) or 
Hp(— k) correspond to massless Dirac Hamiltonians. A HgTe quantum well at d = d c is thus a 
direct solid state realization of a massless Dirac Hamiltonian. Since it does not have any valley 
degeneracy, a d = d c HgTe/(Hg,Cd)Te quantum well is, in a sense, half-graphene. Besides the 
linear term forming the Dirac Hamiltonian, there are additional effects in the HgTe Dirac system, 
such as the quadratic terms in (0Q), and the presence of Zeeman- and inversion asymmetry-induced 
terms, which are discussed in detail in the supplementary material. 

As explained in the introduction, a two-cone Dirac system is the simplest possible realization 
of Dirac fermions for any two dimensional quantum well or thin film, which makes HgTe a very 
interesting model system to investigate Dirac fermion physics. Other benefits include the very 




H D (k) = e(k) + di(k)<7i. 



(1) 



where 



di + id 2 = A[k x — iky) = Ak-, 
d z = M-B{kl + k 2 y ), 
e = C-V(kl + kl). 



(2) 
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high mobility (up to 1.5 x 10 6 cm 2 /Vs for high carrier densities) and the possibility to study the 
effects of a finite relativistic mass Ai (with both positive and negative sign). In this paper, we 
describe magneto transport experiments on gated zero gap HgTe wells that clearly demonstrate 
the Dirac fermion physics expected from Eqs. (0Q) and (O. 

III. EXPERIMENTAL 

For these studies, we have grown by molecular beam epitaxy a number of modulation-doped 
HgTe/Hg 3 Cd .7Te quantum well structures on lattice-matched (Cd,Zn)Te substrates, with a nom- 
inal well width ranging from 5.0 to 7.5 nm (yielding various relativistic masses M.), including 
several samples aiming for the critical thickness of 6.3 nm. From Fig. ED (a), the reader can in- 
fer that the series includes both normal and inverted band gap structures. From X-ray reflectivity 
measurements on our quantum well structures 11 we infer the existence of thickness fluctuations 
of the order of a monolayer in the samples, which corresponds to fluctuations in Ai of around 1 
meV. Subsequently, the wafers have been processed into Hall bar devices with dimensions (length 
L x width W) of (600 x 200) and (20.0 x 13.3) \i m 2 using a low temperature positive optical 
lithography process. For gating purposes a 100 nm thick Si3N 4 /Si02 multilayer gate insulator 
and a 5/50 nm Ti/Au gate electrode are deposited. Ohmic contacts are made by thermal indium 
bonding. A micrograph of such a Hall bar device is shown in the inset of Figf2](a). At zero gate 
voltage, the devices are n-type conducting with carrier concentrations around 5 x 10 11 cm -2 and 
mobilities of several 10 5 cm 2 V -1 s -1 . 

Transport measurements are carried out in a variable temperature magneto-cryostat at a tem- 
perature of 4.2 K, unless indicated otherwise. Typically, a bias voltage of up to 10 mV is applied 
between current contacts 1 and 6 (as denoted in the inset of Fig. [2£a) ), resulting in a current I of 
approximately 1 \x A , as determined by measuring the voltage drop across a reference resistor in 
series with the sample. The resulting longitudinal (V xx , contacts 3 and 5) and transverse (V xy , con- 
tacts 2 and 3) voltages are detected simultaneously yielding the longitudinal (p xx = V xx /I x W/L) 
and transverse (p xy = V xy /I) resistivities. 

Applying a gate voltage V G between the top gate and the 2DEG, the electron density (and 
thus the Fermi energy) can be adjusted. As reported previously, the carrier type can be varied 
from n-type conductance for positive Vg to p-type behavior for negative V G . Hysteresis effects 
due to interfacial states 12 restrict the usable range of gate voltages to |Vq| < 4 V. For reasons of 
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comparison, we have adjusted the gate voltage axes in Figs. 2, 3 and 4 such that V G — V Dirac = 
V corresponds to the Dirac point. V Dirac varies from cool-down to cool-down, but typically is of 
order -1.2 V. 

IV. QUANTUM HALL EFFECT AND THE IDENTIFICATION OF ZERO-GAP SAMPLES 

In Fig. (2a) we plot the Hall conductivity a xy = p xy / (p 2 xy + p 2 xx ) at various fixed magnetic 
fields for a sample with d m d c ~ 6.3 nm as a function of the gate voltage. The conductivity axis 
is correct for the trace taken at 1 T, while the traces for higher fields have been offset by a constant 
amount (in this case one conductance quantum) per Tesla, for reasons that will become obvious 
shortly. First, we note that the traces show well developed quantum Hall plateaus, even for fields 
as low as 1 T. At this low field, the spin-derived Hall-plateaus (the conductance plateaus at an 

2 

even integer times %-) are still less broad than the orbital-induced ones (plateaus at an odd integer 

2 

times %-), which facilitates their assignment. Obviously, because of the large g-factor of HgTe 
(g* = 55.5 for this well, see below) the Landau levels are always spin-resolved. A full assessment 
of Dirac behavior will thus have to come from the field and energy dependence of the Landau level 
structure, which we will provide below. 

First, we will address another question - is the sample really zero gap? Since MBE growth 
calibration is not sufficiently precise to consistently grow a quantum well of exact critical thick- 
ness, we require another independent means to assess the well thickness. We have found a simple 
procedure by analyzing the quantum Hall data of our samples. Specifically, it turns out that the 
crossing point of the lowest Landau levels for the electron and heavy-hole sub-bands is a precise 
measure of well thickness. By solving the Landau levels of the effective Hamiltonian © in a 
magnetic field, we find that each of the spin blocks exhibits a 'zero mode' (n=0 Landau level), 
which is one of the important differences between Landau levels of materials described by a Dirac 
Hamiltonian and those of more traditional metals 13 . The energy of the zero mode is given by 

B$ = C + M-^{V + B), 

E^ = C-M + ^(-V + B) (3) 

for the spin-up and spin-down block, respectively. Here Bj_ is the perpendicular magnetic field. 
The spin splitting, given by Eq — Eq = 2M. — 2£>^-= L , thus increases linearly with magnetic field. 
From ©, we find that there is a critical magnetic field B C L = where the two zero mode spin 



6 



levels become degenerate, Eq = Eq. In the inverted regime Ai/B > 0, this degeneracy occurs 
at a positive magnetic field B\ > 0, while in the normal regime where Ai/B < 0, the crossing 
extrapolates to a negative value of B\ < 0. For a well exactly at the critical thickness d c we have 
Ai/B = 0, and the crossing point will occur at zero field, B\ = 0. Therefore the position of the 
crossing point of the spin states of the lowest electron and hole Landau levels at zero magnetic 
field will give us a direct indication for the existence of a Dirac point (Ai = 0) in the quantum 
well. 

Applying this procedure to the experimental data of Fig. [21(a) (and similar data from the other 
quantum wells in our growth series) is straightforward. Since the Landau levels are already well 
defined at small magnetic fields, we can easily identify the ±1 Landau levels corresponding to the 
two spin blocks of the zero mode as the boundaries of the a xy = plateau, at various magnetic 
fields. The constant offset between the different plots in Fig. [2] (a) implies that we can now 
translate the vertical axis into a field axis with a spacing of 1 T between the scans (the "B"-axis 
in the figure) and we can directly plot the linear spin splitting predicted by Eq. © in Fig. |2](b). 
Extrapolating the linear behavior in the graph allows us to determine B c ± , which in this case leads 
to B\ «0 T- this sample has a Dirac mass close to zero. 

As an illustration of the efficiency and sensitivity of this procedure, Fig. [2] (b) shows the ex- 
traction of Bj_ for three different samples. The sample in the upper panel has an inverted band 
structure since B C L > (from a more detailed fit we find d = 7.0 nm). The middle panel corre- 
sponds to the data of Fig. [2] (a), where the intersection is at B\ ps T, corresponding to Ai = 0, 
and finally the sample in the bottom panel has a not-inverted band structure since the crossing 
point occurs for B\ < 0T (and corresponds to a well- width of approximately 5.7 nm). 

V. FURTHER CHARACTERIZATION OF A ZERO-GAP SAMPLE 

In the following, the sample with B\ = T of Figj2](b) is further investigated. Figs. 12(c) 
and (d) show the Hall conductivity of this sample at 1 and 5 T, respectively, in combination with 
the Shubnikov-de Haas oscillations in the longitudinal resistance. The first thing to note is the 
quantization of the Hall plateaus. Orbital quantization yields plateaus at odd multiples of e 2 /h, 
with additional even-integer plateaus due to spin splitting already observable at 1 T. This is the 
unusual ordering of the Hall plateaus that results from the Dirac Hamiltonian 3 . Moreover, the 
observed plateaus occur at one half the conductance of the plateaus observed for graphen o 14 ' 15 - a 
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direct consequence of the fact that the HgTe quantum well only has a single (spin degenerate) Dirac 
cone, where graphene has two. Furthermore, we always observe a plateau at zero conductance in 
the Hall traces, which is different from the low-field behavior in graphene 3 . The zero conductance 
Hall plateau is always accompanied by a quite large longitudinal resistivity, which is once more 
an indication that - already at IT - the sample is gapped due to spin splitting. 

To further validate our claim that this sample boasts a zero gap Dirac Hamiltonian at low 
energies, we plot in FigfJ] (a) a Landau level fan chart. This chart was obtained by plotting the 
derivative da xy /dY G in a color-coded 3 -dimensional graph as a function of both V G and B±. 
When the sample exhibits a quantum Hall plateau, the Hall conductance obviously is constant 
and its derivative is zero; when a Landau level crosses the Fermi energy, da xy /dWQ reaches a 
maximum, which can be conveniently indicated by the color coding. To translate the gate voltage 
axis to an energy scale for the band structure, we assume that the gate acts as a plane capacitor 
plate, and calculate the electron density in the quantum well as a function of energy using our 
8-band k ■ p model 16 , assuming the well has the critical thickness d = 6.3 nm. Furthermore, in the 
supplementary material, we calculate the density of states as a function of magnetic field for fixed 
electron density and compare the results with the experimental data on the Shubnikov-de Haas 
oscillations. The good agreement of the node position and spin splitting between the experiment 
and theory verifies the validity of the 8-band k ■ p model. The dashed white lines in FigfJ] (a) 
give the Landau level dispersion predicted by our calculation; the very good agreement with the 
experimental peaks in da xy /dWQ is evidence that our Vg to E conversion is self-consistent. 

The Landau-level dispersion in Fig. [3] (a) shows all the characteristics expected from our Dirac 
Hamiltonian (0Q). Besides the zero mode of Eq.©, solving the Landau levels of the effective 

Hamiltonian (OQ) in a magnetic field, yields for the higher Landau levels (n = 1, 2, ) (C, M. = 

0): 

El(n) = ~B ± (2Vn + B) 

+a^2nA 2 ^B ± + (^B^ (V + 2Bn) 2 

Ei(n) = ~B ± (2Dn-B) 
n 

+aJ2nA 2 ^B ± + (~#±) 2 (V - 2Bn) 2 (4) 

where E\, E^ refer to the two spin blocks of our Dirac Hamiltonian, Eq. ©, and a = +, — denote 
the conduction and valence band, respectively. With optimized parameters ( we use C, M. = 
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meV, V = —682 meV-nm 2 , B = —857 meV-nm 2 , A = 373 meV-nm ) the Landau level dispersion 
described by Eq. @i are plotted (dashed white lines) as a function of magnetic field in Fig. [3](b). 
Clearly, the Dirac model agrees well with our experiment for low magnetic field and low-index 
Landau levels, but gradually breaks down when the magnetic field is increased. 

In the low magnetic field limit, one easily finds that Eq. © (for the conduction band) reduces 
to E^\n) w Ay/2n^B± - (~B ± )(2Vn ± B) up to B± linear terms. This corresponds to the 
square-root magnetic field dependence that has meanwhile become the signature of Dirac fermion 
behavior in graphene 3 , with an additional linear term reflecting the large effective g-factor g* 
of the HgTe quantum well. Defining fiB9*B± = £+(0) - £+(0), we find g* w 55.5. There 
are two physical origins for the large g*. Due to the zero gap nature of the present system, the 
most important contribution comes from orbital effects which are fully incorporated in the Dirac 
Hamiltonian. However, there is also a contribution from Zeeman-type terms, which is not included 
in the Dirac Hamiltonian (OQ). This term is less important than the orbital part and will be discussed 
in the supplementary material. 

Another effect that is not included in our model Dirac Hamiltonian is the inversion asymme- 
try of the system. In principle, the HgTe quantum well has structural (SIA) and bulk inversion 
asymmetry (BIA)^, both of which can couple Dirac cones with opposite spin. From the node 
position of Shubnikov-de Haas oscillations (the data are presented in the supplementary material), 
we find that the spin splitting due to SIA is less than 2.5 meV at the largest experimentally acces- 
sible Fermi energy, decaying rapidly with density^. The present experiment does not show any 
evidence of the BIA term; a previous theoretical estimate shows that the BIA term has an energy 
scale of about 1.6 meV 6 . We conclude that also the SIA and BIA terms are small compared to the 
other terms in the Dirac Hamiltonian of Eq. ©. Moreover, they cannot cause the opening of a gap 
in the quantum well spectrum. The relevance of SIA and BIA terms is discussed in more detail in 
the supplementary material. 

VI. ZERO FIELD BEHAVIOR 

Having thus established that we indeed can describe our quantum well as a zero gap Dirac 
system, we now turn to its characteristics at zero magnetic field. Fig. |4](a) plots the resistivity p xx 
vs. gate voltage, often called the Dirac-peak in the graphene community, in this limit. The graph 
clearly shows the expected peaked resistivity and exhibits an asymmetry between n- and p-regime 
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which can be attributed to the large hole mass (increased density of states). In graphene, the width 
of the Dirac-peak is often regarded as a measure of the quality of the sample 18 . The width of the 
Dirac peak in Fig. |4](a) corresponds to a carrier depletion of about An w 3.0 x 10 10 cm~ 2 , which 
is comparable with the situation found in suspended graphene. 

At the Dirac point, we find a minimum conductivity of <J xx ,min = 0.36 e 2 /h at 4.2 K. Its 
temperature dependence is shown in Fig. H](b), which conveys an initially quadratic temperature 
dependence, that for temperatures above about 12 K turns linear. The existence of a finite minimal 
conductivity at vanishing carrier density is a topological (Berry phase) manifestation of the conical 
singularity of the Dirac bands at k = 0. Therefore, our observation of a minimal conductivity in 
HgTe quantum wells provides independent evidence for the Dirac fermion behavior in this ma- 
terial. The observed minimal conductivity (close to e 2 /nhr^-) and the crossover from quadratic 
(oc T 2 ) to linear (oc T) increase with temperature can be understood from calculations based on 
the Kubo formula, in which the current-current correlation function is evaluated for the effective 
Dirac Hamiltonian of Eq. (0Q), assuming the presence of both well width fluctuations and potential 
disorder and including only the dominant terms linear in k. The details of these calculations are 
described in the supplementary material. Qualitatively, the temperature dependence of a xx is 



2e 2 1 
nJl+ (M 2 )/T 2 
e 2 k 2 T 2 



+O^^J, fc s T«r, (5) 
e 2 k T 

&xx oc -r--^-, k B T > T, (6) 
h 1 

where Y is the spectral broadening induced by spin-independent potential disorder. In Eq. © 
the factor of 2 accounts for the spin degeneracy and (Ai 2 ) oc ((d — d c ) 2 ) is the variance of the gap 
due to spatial deviations of the thickness d from the critical value d c . From the X-ray reflectvity 
data on our samples 11 , we estimate a/ (M 2 ) ~ 1 meV. This is comparable with typical values of T 
estimated from the self-consistent Born approximation (see supplementary material). Therefore, 
for a/ (Ai 2 ) ~ T the zero-temperature conductivity is approximately a xx (T — > 0) ~ e 2 /irh, in 
agreement with the data. The T 2 correction reflects the spectral smearing at energies below T. In 
contrast, at k B T > T the linear T-dependence of a xx © reflects the linear density of states, which 
is another manifestation of the Dirac fermion physics in HgTe quantum wells. 

In conclusion, our paper reports the first experimental discovery of a two dimensional massless 
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Dirac fermion in a single valley system. The high mobility in the HgTe quantum wells should 
allow us to directly study ballistic transport phenomena that so far have been hard to access for 
Dirac fermions 21 . Moreover, the material offers an additional parameter for the experiments in 
that the effects of a finite Dirac mass can now be studied in detail. 
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VII. SUPPLEMENTARY ONLINE MATERIAL 

In the supplementary material , we will give theoretical details connected with the analysis of 
the effective g factor, Shubnikov de Haas oscillations and minimum conductivity in HgTe quantum 
wells. 

A. Effective Hamiltonian for HgTe quantum wells 

In this section, we will discuss the complete expression of the effective Hamiltonian of HgTe 
quantum wells near the critical thickness d c . As first described by Bernevig, Hughes and Zhang 4 , 
the low energy physics of HgTe quantum wells is determined by four states \Ei, |), \H ly |), 
\Ei, — |) and \Hi, — §). With these four states as basis, the complete Hamiltonian of the system 
when a magnetic field B± is applied in the the z-direction can be written as 

H = H e ff + Hzeeman + H SIA + H BIA (7) 

The effective Hamiltonian is given by Eq. (1) in the main part of the article, where a Peierls 
substitution k — > k + | A has been applied (A is the magnetic vector potential). All parameters 
A, B, C, V and At can be determined by fitting to the experimental Landau level dispersion at 
T = 4.2 K for u^deg (we neglect in our calculations the dependence of the parameters on 
the electron density), which are listed in the table HI 
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The Zeeman term H Ze 



has the form 6 



/ 



H 



Zee 



9e 



\ 



9h 



-9e 



\ 



(8) 



-9h J 



with hb — 2^ an d ^ f° ur band effective g-factor g& for E\ (Hi) bands. 

Since in HgTe quantum wells the inversion symmetry is broken, we also need to discuss two 
terms that result from inversion asymmetry, i.e. the structural inversion asymmetry (SIA) term and 
the bulk inversion asymmetry (BIA) term. The SIA term is due to the asymmetry of the quantum 
well potential and has the form 17 
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(9) 
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In (001) grown HgTe quantum wells,, the Rashba spin splitting is dominantly responsible for 
the beating pattern of the Shubnikov-de Haas oscillations. By comparing the results of a Kane 
model calculation with the experimental data, one can thus directly determine the Rashba coeffi- 
cient. As discussed in the following section, we find that for a quite large range of gate voltages, 
the Rashba spin splitting is less than 2.5 meV near the Fermi energy, which corresponds to £ e ~ 16 
meV-nm, x ~ 2.0 meV-nm 2 and £ h w 5.0 meV-nm 3 . Furthermore we find that the electron Rashba 
splitting £ e is always dominant over the other two terms. 

Since the zinc-blend crystal structure of HgTe is not inversion-symmetric, additional BIA terms 
appear in the effective Hamiltonian, given by 6 

/ -A \ 
A 
A 
\-A / 

Since the BIA Hamiltonian is a constant term, it will only change the the Dirac point to a circle and 
the system remains gapless. An early estimate of magnitude of BIA term gives A « 1.6meV 6 , 



H 



BIA 



(10) 
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which is of the same order of magnitude as the disorder broadening and the fluctuations of the band 
gap due to variations in well width. Therefore in the present experiment we do not find evidence 
of a pronounced effect from this term. 

TABLE I: The parameters used in the effective model. We take C to be zero to shift 
the Dirac point to zero energy. 



„4.(meV • nm) 


373 


£>(meV • nm) 


-857 


P(meV) 


-682 


M(meV) 


-0.035 


A (meV) 


1.6 


9E 


18.5 


9H 


2.4 



Neglecting the SIA and BIA terms, the Landau level spectrum is described by 

El{n) = - G -^(2Vn + B) + ^(g E + g H ) 

+a ^2nA^ +{M-B ± (|0D + 2Bn) - ^f{g E - g H ))) 2 

El(n) = -^Vn - B) - ^{g E + g H ) 

+a^2nA^ + (m - B ± + 2Bn) + ^f(g E - ^))) 2 , 

(11) 

where n = 1,2,..., and a = + (a = — ) for the conduction (valence band) and the parameter C 
is taken to be zero, setting the Dirac point at zero energy. The zero mode states (n — 0) have the 
dispersion: 



El = M- e -^{V + B) + ^ 9E 

Ei = -M + ^(-V + B)-^g H , (12) 



and the zero mode splitting is: 



AE s = 2M- 2 ^B + ^(g E + 9H ). (13) 
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Finally, we note that the total spin splitting, defined as the energy difference between the spin- 
up and spin-down zero modes, has two origins in the four band effective model (A4 — >■ 0). 
One comes from the Zeeman term, which gives the energy splitting AE s2 = //b^-^J- = 
P B 2 x (g E + g H ) with g\ = 10.5. The other origin stems from the combined orbital effects of the 
linear and quadratic terms in the Hamitonian H e ff, and is given by AE s i = [iBB±gl = —^j^-B 
with gl = 45. These two terms together give the effective g-factor g* defined in the main part of 
the article. 



B. Shubnikov-de Haas oscillations 



In order to compare the Kane-model calculations with the experimental data, the density of 
states (DOS) at the Fermi level was calculated from the Landau level spectrum (see Ref. 16 for 
more details). The Shubnikov-de Haas (SdH) oscillations observed in the experiments are directly 
related to the oscillations of the DOS at the Fermi energy. In Fig.[5]the calculated DOS, broadened 
by convolution with a Gaussian with a width r =1.2 meV, is displayed together with SdH data 
for three different values of the gate voltage. A very good agreement between experiment and 
theory is evident. For all three different gate voltages, we find that there are always two sets of 
minima in the oscillations, one deep and one shallow, which result from the two sets of Landau 
levels for opposite spin (cf. Eqs. (fTTI) and (PT21)). We first make the obvious identification that 
the deep minima result from the Landau level splitting and the shallow minima from the spin 
splitting. It is now instructive to plot the inverse of the magnetic field value for the positions of 
the deep minima as a function of the number N associated with the deep minima of the oscillation, 
as done previously to demonstrate the implications of the Berry phase of the Dirac Hamiltonian 
for the quantum Hall effect in graphene^. As shown in Fig. 6, we find that when V G = 2 V, 
a straight line fit to the data extrapolates to N = 0, while for Vg = V, the fit extrapolates to 
iV = 1/2. This different behavior can be understood from the effective Hamiltonian (1) in the 
main text of the article. For Vq = V, the electron density is low and the Fermi energy is near 
the Dirac point. In this limit, the band dispersion is dominated by the linear term in wave vector 
k and the deep minima correspond to the filling factors u — 2(N — |) (the factor of 2 takes the 
spin into account). The intercept of the straight-line fit evidently corresponds to the filling factor 
v = 0, which explains the N — 1/2 intercept. However, for Vg = 2 V, the electron density is 
increased and the Fermi energy is far away from the Dirac point. Consequently, the linear term is 
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no longer dominant and other terms, such as the quadratic ones, will come into play. In this limit, 
the system recovers the usual behavior of a two dimensional electron gas and the deep minima 
correspond to the filling factors v = 2N, implying that the intercept occurs at N = 0. Note 
that in HgTe, in contrast with graphene, one always has additional shallow minima besides the 
deep ones due to the coexistence of the linear term and other type of terms, such as quadratic or 
Zeeman terms. Another important feature of the SdH oscillations in Fig. [5] is the appearance of a 
beating pattern when V G = 2V, indicating the occurrence of Rashba spin splitting 22 , which from 
this data is estimated to be 2.5 meV. The beating feature is not observable when the gate voltage is 
in the range of ~ 1 V, which indicates that the system becomes more symmetric for low electron 
densities. An extensive discussion of this effect can be found in Ref. 17 . 



C. Calculation of the minimal conductivity 

In this section of the supplementary material, we discuss details of the calculation of the mini- 
mal conductivity given by Eqs. (5) and (6) in the main part of the article. 
We use the Kubo formula for the longitudinal (xx) dc conductivity, 



df\ r d 2 k 



v x A k , v x A k . 



(14) 



U A.R 



ik( ^- G M (15) 

ft F)h ' 9-7T7 

Here the velocity operator v x , spectral function A k e and retarded/advanced Green's functions 
are 2x2 matrices in Ei-Hi subband space of the effective Hamiltonian H o(k) described by 
Eqs. (1) and (2) of the main manuscript. We use the symbol (Tr) to designate the trace operation 
. The spin degree of freedom is accounted for by the factor of 2 in Eq. (fl4l) . /(e) is the Fermi 
function, e is the energy measured from the neutrality point and k = (k xi k yi 0) is the wave vector 
in the plane of the quantum well (QW). 

Our next step is to calculate the disorder-averaged Green's functions G ke . In the minimal 
conductivity regime the most relevant types of disorder in our HgTe QWs are the inhomogeneity 
of the carrier density and the spatial fluctuations of the QW thickness d(r) around critical value 
d c . The carrier density inhomogeneity induces random fluctuations of the electrostatic potential in 
the QW, which we treat as weak gaussian disorder with standard averaging procedures leading to 



15 



the complex self-energy in the equation for G^f (see, Eqs. (fT9l ) and (1201) below). In this respect, 
we follow the previous theoretical work on graphene (see, e.g. Refs. 2 - <^). The new feature of our 
model is that it also accounts for the QW thickness fluctuations which is a specific type of disorder 
in HgTe/CdTe structures. This type of disorder induces sizable regions in the sample where the 
effective Dirac mass Ai has small positive or negative values oc d(r) — d c : 

M(v)^M'-{d(r)-d c ), (16) 

where Ai' ss 15 meV-nm -1 is a proportionality coefficient that we determine from band structure 
calculations. We assume that the gap Ai(r) varies slowly between the regions with Ai > and 
Ai < in the sense that the carrier motion adiabatically adjusts to the variation of Ai(r). Both 
Ai (r) and its gradient V r Ai (r) are supposed to vanish upon averaging over the whole sample 
area (a) so that the leading nonzero moment of Ai (r) is the variance: 

(M{t)) = (V r M(r)) = 0, (M 2 (r)>^0, (...) = a' 1 J dr ... (17) 

In view of the adiabatic dependence Ai (r) it is convenient to use the mixed representation for the 
Green's functions defined by the Wigner transformation: 

O) = / «3f* (r + y. ' - y) e- ik - (18) 

where G^(r) satisfies the equation: 

Here we omit the k 2 corrections to the linear Dirac Hamiltonian [see Eqs. (1) and (2) for H D (k) 
in the main manuscript] because the main contribution to the minimal conductivity comes from 
the vicinity of the k = point (fj, is the Fermi energy measured from the neutrality point). For the 
same reason, in Eq. (fT9l ) the self-energy S^ /4 (r) (generated by the random potential fluctuations) 
is taken at k = 0. It has been established earlier that the universal minimal conductivity follows 
already from the self-consistent Born approximation or equivalent approaches (e.g. Refs.— ). 
We also adopt this approximation for the self-energy: 

E <>> = /||C.Or). C - q = C<„ (20) 



XL r 



Act k 



V r - a z Ai 



:Vi 
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where £ q is the Fourier transform of the correlation function of the random potential, which is an 
even function of the wave-vector q due to the statistical homogeneity of the disorder. 

In order to solve Eq. (fl9l) we follow the same strategy as in the case of the uniform Ai, i.e. 
we first apply operator e + fi — Sf M (r) + Act (k — \ V r ) + o z M(r + §Vk) to both sides of the 
equation from the left. Since k — |V r does not commute with M(r + f Vk) and Ef /A (r), there 
appear additional gradient terms which are absent if M.(v) = const: 



e + fi - Xr>)) 2 -M 2 (r + |V k ) - A 2 (k - |V r )' 
+^<r (zV r Ef A (r) + z x V r M(r) + z x V r .M (r + |V k ) 



(21) 



e + fi 



r 



Acrk + a z M(r). 



We are interested in the average Green's function over the sample area: G^ /A = a 1 J G^ A (r) dr 



and, respectively, E 



,-1 



J E^ M (r) dr. Upon such averaging the linear terms Ai(r), 



V r .M(r) and V r E^ /A (r) in Eq. (I2TI) vanish, while the quadratic term M 2 (r + fVk) does not. 
The latter is assumed to vary slowly in space, allowing us to neglect the corrections oc Vk and to 
obtain the following equations for the Green's function and the self-energy: 



^r/a _ e + /i - Ef M + Acrk R/A _ f d 2 q( q e + n - Ef /A 

k ' e ~ (e + ^u-Ef- 4 ) 2 - (M 2 ) -AW e " J (2tt) 2 (e + /i - Ef M ) 2 - (M 2 ) - A 2 q 2 ' 

(22) 

These results are justified if the average slope of the gap variation, a/ (| V r M | 2 ) / a/ {M 2 }, is small 
compared to the characteristic electron wave-number T e /A, 



W r M\ 2 ) T e R/ A 

< -r, T £ = |ImE e |, (23) 



where T e is the spectral broadening due to the finite elastic life-time. It eliminates the infrared 
divergence of the k integral in Kubo formula (fl4)) . providing an effective cut-off at small values of 
k. 

We now use Eq. (1221) for G^ /A and v x = dH D (k)/hdk x « (A/H)a x to calculate the k integral 
in Eq. (fT4l) and express the conductivity in the form: 

r°° de 

°~ {T) = L w^ps a - w " (24) 
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(TazrAE) 



2e^ 1 
Vh X 2 



e 2 + T 



2|e|r f 



arcsm 



+ 



2|e|r e " v/( r e-e 2 + (-^ 2 )) 2 + 4e 2 r2 

(r 2 + e 2 ) 2 + (r 2 -e 2 )(.M 2 >l _ 

e = e + /i — Re L 



(25) 



(r 2 -e 2 + (M 2 )) 2 + 4e 2 r 2 
For e, /i — > Eq. (|25l) yields the zero-temperature minimal conductivity a xx (0). To completely 
specify Eq. (|25l) we find the self-energy from Eq. (1221) in the form of a power-law expansion: 



Co + C ie + C 2 e 2 + ... 



(26) 



Assuming Cq = Co = const and a cutoff A = A 2 /B ~ 120 meV at high energies [where 
the quadratic term Bk 2 becomes comparable with the linear one Ak in ifo(k)], we obtain the 
expansion coefficients as 26 



C = TtT, r = v/A 2 e- 2 /« - (M 2 ), a = Co/2nA 2 , (27) 
Ci = - (r^) Tf)'\c 2 = ±» a ff + -%ffl Cg, u/r«a«l. (28) 

This model adequately describes the observed minimal conductivity and its temperature depen- 
dence. In particular, at low temperatures ksT <ti Y the T-dependence is quadratic: 

1 



trr* 2e 2 

&xx{T) w — 



l/9+(M 2 )/r 2 n 2 k 2 B T 2 
~\~ — , _ — „ x 



(29) 



tt/i [i + (M 2 )/r 2 (i + (M 2 )/r 2 ) 3 r 2 

At ksT > Y the conductivity becomes approximately linear as a manifestation of the linear density 
of states of the 2D Dirac fermions in HgTe QWs. The fit to the experimental curve a xx {T) in 
Fig. 4b is achieved for a/(M 2 ) = Y = 1 meV and fi/a 3 Y = 0.06. 
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FIG. 1: (a) Sub-band energies of HgTeZHgo.3Cdo.7Te quantum wells as a function of well thickness d. (b) 
In plane dispersion of a quantum well at the critical thickness d c ~ 6.3 nm. (c) A 3D plot of the Dirac cone 
describing the low energy spectrum for d = d c . 
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FIG. 2: (a) Gate voltage dependence of the transverse conductivity for a HgTe quantum well with a well 
thickness d w 6.3 nm for fields of (from top to bottom) 1 (purple), 2 (green), 3 (blue), 4 (red), and 5 (black) 
Tesla. Note that the conductivity axis belongs to the 1 Tesla trace. All traces for higher magnetic fields 
have been shifted down by one conductance quantum per Tesla so as to show the determination of Bj_, 
using the "B"-axis on the right-hand ordinate as field axis (as based on Eq. [3l see text). Inset: Micrograph 
of a typical sample, (b) Determination of for three exemplary quantum wells grown close to d c . By 
extrapolating the crossing of the zero level spin states to the gate voltage needed to reach the Dirac point 
one finds B C L ; when B c ^ = the well has the critical thickness d c ss 6.3 nm. See text for details. (c,d) The 
transverse conductivity and longitudinal resistivity of a sample with well thickness d m 6.3 nm at 1 and 5 
Tesla, respectively. 
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FIG. 3: Experimental Landau level fan charts obtained by plotting da xy /dNQ in a color-coded 3- 
dimensional graph as a function of both Nq and B±. Energies are measured with respect to the Dirac 
point. In (a) this data is compared with the results of a calculation (dashed black-and-white line) from our 
8-band k-p model, (b) Comparison with the fan chart computed from the Dirac Hamiltonian (Eq. [T]) (dashed 
black-and-white line) using the parameters described in the text. 




FIG. 4: Resistivity p xx vs. gate voltage Vq (Dirac -peak) measured on a zero gap quantum well in the ab- 
sence of an external magnetic field. This trace was taken at a sample temperature of 4.2 K. (b) Conductivity 
at the Dirac point as a function of temperature. Open circles are experimental data, the red line is a fit to 
Eqs. (5) and (6). 




FIG. 5: Calculated Fermi level density of states of a d = d c HgTe quantum well for various positive gate 
voltages as a function of perpendicular magnetic field (thin red lines) compared with the experimental SdH 
oscillations (black thick lines). 
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FIG. 6: The inverse of the magnetic field at which the deep minima in the Shubnikov-de Haas oscillations 
occur, as a function of ordinal number N for Vq = 2 V (red circle) and Vg = V (blue circle). The ordinal 
number N can also be related to the conductivity at the corresponding Hall plateau. For Vq = V, the 
corresponding Hall conductivity for deep minimum N is given by a xy = ^-2(N — |), while for Vq = 2 V, 

2 

the Hall conductivity is a xy = jr2N. 



